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Abstract 

We use the two dimensional extended Hubbard Hamiltonian with the position 
of the attractive potential as a variable parameter with a BCS type approach 
to study the interplay between the superconductor transition temperature T c 
and hole content for high temperature superconductors. This novel method 
gives some insight on the range and intensity of the Cooper pair interac- 
tion and why different compounds have different values for their measured 
coherence lengths and it describes well the experimental results of the super- 
conducting phase diagram T c x n. The calculations may also be used to study 
the effect of the applied pressure with the assumption that it increases the at- 
tractive potential which is accompanied by an increase in the superconductor 
gap. In this way we obtain a microscopic interpretation for the intrinsic term 
and a general expansion for T c in terms of the pressure which reproduces well 

the experimental measurements on the S^S^CaC^Os+ySuperconductors. 
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I. INTRODUCTION 



In order to understand the superconductivity mechanism for high-T c superconductors 
(HTSC), many experiments have attempted to find out possibles correlations between the 
superconductivity temperatures T c and some physical parameter. Some of this experiments 
deal with the interplay among T c , hole contents and pressure effects which have revealed 
interesting features like the maximum attained T c M and have attracted a lot of attention as 
it is summarized by some review articles On the other hand, neither one of the several 

microscopic theoretical models have been able to deal successfully with all these observed 
basic properties and therefore the nature of the condensate pair has not yet been established. 
As a consequence, some phenomenological approaches have been advanced in an effort to 
account for the large experimental data on HTSC. For instance, the T c dependence on the 
number of carriers per CuO^ follows the so-called universal inverse parabolic curve, i.e., 
T c = T™ ax [l — r]{n — n op ) 2 } where n op is the optimum carrier density n and i] is an adjustable 
parameter which depends on the type of compound. 

Recently we have introduced an intermediate approach |5j which is neither pure mi- 
croscopic nor pure phenomenological: it starts with the microscopic extended Hubbard 
Hamiltonian (t-U-V) on a square lattice of side a 

H = - *( c L c j«r + h - c -) + U J2 n H n iT ~ V H n i n r W 

where t is the nearest neighbor hopping integral, U is the on-site repulsion and V a phe- 
nomenological attractive interaction. Several authors have argued that this Hamiltonian 
accounts for the physics of the charge carriers |||7j in HTSC. The strength of the phe- 
nomenological potential coupling constants are determined a posteriori in comparison with 
the experimental results for T c x n. The hopping integral t is obtained by comparing with 
band structure calculations [0. Using the BCS method and the Fourier transform of the po- 
tential we have derived zero and finite temperature gap equations which are integrated up to 
region where pair formation has a non-negligible probability. In this way we have been able 
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to relate T c , hole content n and the strength of the zero temperature gap A(0) as described 
in detail in Ref. 0. As concerns the phase diagram, its optimal doping is very sensitive if 
we move V from its original nearest neighbor position for the s x 2 +y 2 -channel but does not 
change its optimal position for d x 2_ y 2-wave. A study of these and others channels as well 
as their complex mixtures and interplay as function of the position of V will be published 
elsewhere ||. This method was successfully applied to the experimental measurements of 
the curves T c xn for the La(Sr, Ba)CuO and Y BaCuO compounds f|. On the other hand if 
we study the typical parameter which characterizes the strength of the attractive interaction 
for a square well in 2D with range d and strength §| V, x = m e d 2 V/2h 2 , we see that for the 
above two compounds and Hgl201 and Bi 2 Sr 2 CaCu20s+y, it varies between 10-3000 and 
such high values indicate that the size of the bound states (or Cooper pairs) is related to the 
position of the minimum and the range of the attractive potential. This is the main motiva- 
tion to change the position of the attractive potential V from the usual nearest neighbor and 
it is very interesting to find out that the curves for T c xn yield a different s^+^-wave optimal 
doping value depending on the size of V while it remains fixed for d x 2_ y 2-w&ve. Another in- 
teresting point from the theoretical point of view is that it provides a natural interpretation 



for the order parameter expansion introduced by the spin-fluctuation theory [IU|. In this 
approach it has been proposed a d x 2_ y 2-w&ve gap of the form Ao ( k) (cosk x a — cosk y a) where 
A (k) is expanded in powers of (cosk x a + cosk y a). Analising this expansion, one can easily 
verify that it contains terms like (cosk x a — cosk y a)(cosk x a + cosk y a) which are proportional 
to (cos(2k x a) — cos(2k y a)) and what can be seen as a type of d-wave gap symmetry that 
arises from a potential like V(k) = Vo(cos(2k x a) + cos(2k y a)). By the same token, we can 
find terms proportional to (cos(3k x a) — cos(3k y a)) which can be originated by a potential 
V(k) = Vo(cos(3k x a) + cos(3k y a)) and so on. Based on this result, we adopt below a notation 
V(k) = V (cos(bk x a) + cos(bk y a)), with b = 1,2,3, .... As one expects, the more distant is 
the position of V, that is the larger is b, smaller is the value of the density necessary to 
start the superconductivity process of pair formation. In this way we have derived T c x n 
curves H that agree well with the experimental data when V was placed at 6 th (s-wave) and 



3 rd (d-wave) neighbors for La(Sr, Ba)CuO and YBaCuO compounds respectively which by 
the way is very close the ratios of their measured coherence length £ [[II] . Thus the method 
provides a very interesting and direct interpretation for the measured values of the coherence 
length for the HTSC compounds besides the fact that it gives a hint for the superconducting 
fundamental interaction in HTSC. 



We have also extended this approach to study the pressure effects |L2[ on HTSC. We have 
noticed that T c (n) is very sensitive to changes in the zero temperature gap amplitude A(0) 
which is used as an adjustable parameter. Consequently we have introduced the assumption 
that the pressure induces changes on the attractive potential end ultimately on A(0) ||12|| . 
This assumption together with the well documented fact that the applied pressure induces 
a charge transfer from the reservoir layers to the Cu0 2 planes, has enabled us to derive an 
expansion for T c (n, P) in power of the pressure. This expansion was applied to the Hgl201, 
1212 and 1232 compounds and has reproduced the measured values of T c for a large set 
of hole content and applied pressures up to 45 GPa |I2]| . We should stress that since the 
principal interaction mechanism which forms Cooper pairs has yet to be identified, the 
connection between T c , hole content, pressure effects and a microscopic Hamiltonian should 
provide some clues for this interaction and it is a efficient method to obtain information on 
the microscopic mechanism for HTSC. Thus our study provides a guide to the strength and 
range of this fundamental mechanism. 

In this paper we shall clarify and illustrate the main physical ideas of the method men- 
tioned above and to show that it is general and can be applied to other HTSC than the 
Hg-based compounds, as we apply it to the experimental phase diagram and pressure data 
on Bi2212. The Bi 2 Sr 2 C 'aCu 2 Os+ y family has been studied by several groups and in partic- 
ular, we shall discuss firstly the measurements of the dependence of the critical temperature 
T c on the hole content n [T3] to set the pertinent parameters to be used afterwards on the 



study of the variation of T c under applied external pressures |14]]. The Bi 2 Sr 2 C 'aCu 2 's+yis 
strongly anisotropic, its resistivity on the ab plane is 10 4 larger than its similar along the 



c-axis 



14| which justify the use of a 2D model and has one of the shortest coherence length 



Il| among all cuprates superconductors and this property will be taken into account to set 
the parameter b in our calculations. 



II. THE METHOD 



Following the BCS method [15] we use a many-body wave function which describes the 



formation of pairs of charge carriers 

|$) = n(^ + ^at T a^)|$o) ) (2) 



A- 



where |<3> ) is the empty band state and w| + t>| = 1. 

In connection with the Hamiltonian of Eq. ([!]), and the variational approach we 
obtain 

*s = - E ( 3 ) 

I 

where Vzr is the Fourier transform of the potential of Eq. ([!]), and following 
along the lines of Ref. |7j may be written in a "separable" form, V-g- = U — 
2V (cos(k x a) cos(l x a) + cos(k y a) cos(l y a)). A similar approach which takes next-nearest 
neighbor hopping into account was developed by Angilella et al J7| with the same purpose 
to study pressure effects on Bi2Sr 2 C aC^Os+yCompouiads but they only deal with nearest 
neighbor attractions. With the use of the expressions for Uf and v g, we get an equation for 
the zero temperature gap, 



r 2 + Af 



Where ^ = —2t(cos(k x a) + cos(k v a)) — /i , 

Thus since the gap has the same functional form of the potential, namely, = 
A (0) (cos(k x a) ± cos(k y a))/2, where the plus sign is for the s-wave and the minus sign is 
for d-wave channel. The chemical potential \i yields the density n that must be calculated 
self-consistently |§ but as it concerns the superconducting phase boundary it suffices to 
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approximate it by the T = OK value of the maximum energy (concentration dependent) in 
the tight-binding band. 

In order to study pair formation, we calculate the probability of finding a hole pair, that 
is, the condensation amplitude = uzv^. It is a general result ]nj that it has a maximum 
at hp and drops very rapidly for \k\ > kp. According to these considerations and taking 
k = (0, 0) in Eq. (f|), we obtain for s-wave 

J " Jo ((2i(/ + ( a ,/3)-^)) 2 + (A(0)/ + (6a,6/3)/2)2) 

where /±(a,/3) = cosa ± cos/3 and a = k x a and /3 = k y a. \i = f + (a F ,f3 F ) and ap and 
/3^ are the maximum T = occupied values (like a Fermi momentum) that depend on the 
density of holes n. The parameter b is used to set the position of the attractive potential V, 
thus b = 2 for third nearest neighbor and so on. We use ap — [3p ~ arcos(l — n) and 
/i = — 2t(cos(ap)). The integrations are performed up to oa and /3a, which are determined 
at values where the condensation amplitude becomes very small, namely, w 0.01. This 
is usually attained for £jj > 6A(0). 

For T^O, the excitations with their respective probability must be taken into account. 
The derivation of a self-consistently temperature-dependent gap equation is analogous to 
that which leads to Eq. (|^). At this point, we again follow the BCS approach |15| and assume 



that A(T) vanishes at the critical temperature T c , which yields the following equation for 
s-wave symmetry, 

„ r » ffA (!7-2l// + (f> a ,V))/ + (KV)tanh(J|fj) 

B = l i ^ 2t( Ma ,S)- U ( aF ,M) ' <6) 

where we again integrate up to a& and /3a 

Now, if we had integrated Eqs. (|5] and 0) over the whole Brillouin Zone and solve them 

self-consistently, we would get A = B = — (2tt) 2 . Instead of this procedure, we find, for a 

given density n, the value of T c which is compatible, i.e., make the equality A — B for a 

given A(0) and b. This is the basic procedure of our method. The same approach can be 

carried out for d-wave symmetry but in this case we take k = (0,7r/2) in Eqs. ([5] and [|). 



Thus, with only one given value for the parameter b and another one for A(0), we can find 
the value of T c for a given density of carriers n and therefore obtain the curve T c x n. 



III. PRESSURE EFFECTS 

Now we turn our attention to the study of the HTSC under pressure. One of the effects 
of the pressure which is generally accepted and well documented in certain materials is 
an increase of the carrier concentration on the CuO<i planes transferred from the reservoir 
layers 0-|4j]. Such pressure induced charge transfer (PICT) has been confirmed by Hall 
effect and thermoelectric power measurements on several compounds ||. Therefore this 
effect combined with an assumption of an intrinsic variation of T c (linear on the pressure) 
independently of the charge transfer was largely explored to account for the quantitative 



relation between T c and the pressure P and it gave origin to many models [16-20]. However, 



the various pressure data taken with the mercury family of compounds pT|,|I| for underdoped 
and overdoped compounds of HgBa^CuO^^ (Hgl201) has some features that could not be 
interpreted piyifl by these pure charge transfer models by a single set of parameters. On 
the other hand we could explain the Hg-family data using a single set of parameter and this 
is an achievement of our method described below. 

The study of T c x n which is necessary to find the appropriated values of b and A(0) (at 
P = 0) demonstrated the direct proportionality between A(0) and T c which varies with the 
hole content n and it is maximum at the optimal doping. Taking this result into account 
and the discussion of the above paragraph, we are led to propose that the effects of pressure 
are two-fold: (i)- The well accepted PICT and an additional assumption; (ii)- A change in 
the attractive potential V as proposed previously by Angilella et al |7[ which also implies 
on a change in the zero temperature gap A(0). The assumption "ii" came from the analysis 
of the curves T c x n which have revealed that T c is proportional to A(0) but with a large 
constant of proportionality around the optimal doping [0 n op and smaller near the both 
underdoped and overdoped extremes of the T c x n curve. The quantitative changes on V 



S 



due to the applied pressure were estimated in Ref. [0] for the Bi^Sr-iC aC^Os+j/Using the 
measured values |14| of the compressibility tensor under the linear approximation and the 



self-consistently equations. 

Thus the PICT (i) implies that n(P) = n + An(P) and the assumption of a pressure 
dependent gap (ii) implies that A(0, P) = A(0) + AA(0,P) and both equations lead to 
T c (n, P) = T c (n(P), A(0, P)). Therefore to estimate T c for a given compound with a nominal 
value of n and under a given pressure P, we perform an expansion of T c (n, P) in terms of 
P. With the assumption of the linear variation of n and T™ ax (or A(0)) on the pressure, we 
obtain only terms up to third order, that is, 

T c (n,P)= j2a z P z /Z\ (7a) 
z=o 

with 

az = i dW)^ 1 + ikw )ZTMn m p)) > z = h 2 and 3 (7b) 

and ao == T c (n(0), A(0)) is just the P = term and where the derivatives d/ ^y and |p 

dP 



are determined by comparing with the experimental data; d/ tp can be estimated by a set 



of values T v x P for a given density n and can be determined by Hall coefficients or 
thermoelectric power measurements. After determined these parameters, one can derive 
analytical expressions for each coefficient as function of n; using the universal parabolic 
fitting we can obtain a analytical expression for dT c /dn and studying the phase diagram for 
different values of A(0), we can calculate 8T c /dP which calculates the T c dependency on 
P (assumption ii). This procedure gives an intrinsic term which has a clear interpretation 
since it comes from the changes in A(0) as well as a new third order term. It is remarkable 
that with a single set of parameters calculated for HG1201, with just the change in T™ ax , we 
were also able to apply the same expansion (with same coefficients) to HG1212 and HG1223. 
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IV. COMPARISON WITH THE EXPERIMENTS 



In order to apply the entire approach described above, we start with the P = curve 
T c xn for Bi 2 Sr 2 CaCu20^ +y which allows us to determine the initial parameters b and A(0). 
We have used the values t = 0.05eV and V= 0.052eV given in Ref. 0. The coherence length 
estimate for Bi 2 Sr2CaCu20s+ y is ~ 13A which is about 10% less than the £ for YBaCuO 
fTT| . Since we have used b = 2 — 3 for fitting the phase diagram for these compounds, we use 
here b = 2.6 — 2.7 in order to obtain T c xn for i^SV^CaC^Os+y. The results are displayed 
on Fig.l and one can see that it provides a very reasonable fitting for the experimental 



points of Allgeier et al ||13|| . In order to match the T c max = 92.5 with V= 0.052eV, we used 
A(0) = 214K which is the same for all values of n. This value is very close A(0) = 210K 
used to the Hgl210 compounds because both compounds have very closed values for T c ma:r = 
92 — 94K and we used for both a s x 2 +y 2 gap function. We have also plotted a parametric 
curve which has been applied to fit the experimental points, T c = T c ma:E [l — r](n — n op ) 2 ] 
with T c max = 92. 5K and 77 = 9.0. However the data of Huang et al [[14] appears to have 
smaller T™ ax , displaced toward higher hole densities and they are fitted with rj = 25.0. The 
parabolic is a pure phenomenological fitting which is used here only because it provides a 
analytical expression that can be used in the pressure expansion to estimate dT c /dn. 

Now that b and A(0) has been determined let us study the Bi2Sr2CaCu20$ +y compo\md 
under applied pressure. From the analysis of the phase diagram shown in Fig.l and 
by varying the value of A(0) around the P = value of A(0) = 214K, we can infer that 
dT c /dA(0) = 0.35 at n « 0.25 which is near the optimal doping value of hole content used by 
Angilella et al 0. From the Hall coefficient measurements fl4jl , we get d(ln(n))/dP = +8% 
and since their results seens to be shifted to larger values of hole contents, we estimate 
dn/dP ~ 3.0 x 10 -2 which is two order of magnitude higher than previously derived dn/dP ~ 
1.8 x 10~ 3 for the pressure effects on Eg compounds From the low pressure data we can 
get «i = Q.OK/GPa which allows us to derive the value dA(0)/dP = 4.86K/GPa which 
is necessary to determine all the high order coefficients. We noticed that this value is very 
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close to dA(0)/dP = 4.30K/GPa derived for the Hg-based HTSC. Now , with the values 
of dA(0)/dP and dn/dP determined, we can calculate the high order terms of Eq.|75| and 
in particular, we obtain 02 = 2A8K/GPa 2 . The third order coefficient was also calculated 
but it does not give any appreciable contribution in the range of pressure below 2GPa. 
Thus with all the coefficients of Eq]75| determined, we are all set to calculate the values 
of T c (n,P). The calculated results are in excellent agreement with the experimental data 
14j as it is shown in Fig. 2. The large value of dn/dP yields a very large negative second 



order coefficient «2 what makes the curve change its derivative at P « 1.2 GPa. For the Hg 
compounds the change in sign of dT c /dP occurs at much larger P « 30.0 GPa and from the 
above analysis we can see that this difference in behavior is due to the large difference on 
the value of the charge transfer dn/dP for this two compounds. 

V. CONCLUSIONS 

Thus, we conclude this work pointing out that our novel calculations based on a BCS 
type mean field and on the extended Hubbard Hamiltonian with the position of V as a 
variable parameter, is appropriate to describe the interplay between the T c and hole content 
for HTSC. The change in the position of the attractive interaction provides informations 
on the range of the Cooper pair attractive mechanism. Our method yields also a novel and 
systematic way to study the effects of the pressure and gives indication that the pressure 
induces a variation on the attractive potential which in turn, gives an intuitive interpreta- 
tion on the origin of the intrinsic term. We should also emphasize that our procedure is 
based and has its starting point on a microscopic Hamiltonian and it differs from some pure 
phenomenological ways to fit the data 2!| that uses a parabolic fitting to the experimental 



pressure data and are in the same category of the inverse parabolic fitting for the T c x n 
curves mentioned above that are useful to obtain information on the different sets of data 
but has not any microscopic implications which is our main concern here. By applying it to 
the i^S^CaC^Og+ySuperconductors anc i comparing the results with our previous calcu- 
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lations on Hg-based HTSC we have seen that the method not only deals successfully with 
the the experimental results but also displays clearly the physical properties that accounts 
for their different experimental behavior under pressure. 
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FIGURES 



Fig.l Calculations for the Bi phase diagram. The continuous line is the phenomenological 
fitting and the dashed line are the calculations described in the article. The squares are 
experimental points taken from Ref.13. 

Fig. 2 Calculations for the variations of T c for the Bi compound under pressure. The circles 
are experimental points taken from Ref.14. 
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